Thermally-activated charge reversibility of gallium vacancies in GaAs 
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The dominant charge state for the Ga vacancy in GaAs has been the subject of a long debate, 
with experiments proposing — 1, —2 or —3 as the best answer. We revisit this problem using ab 
initio calculations to compute the effects of temperature on the Gibbs free energy of formation, and 
we find that the thermal dependence of the Fermi level and of the ionization levels lead to a reversal 
of the preferred charge state as the temperature increases. Calculating the concentrations of gallium 
vacancies based on these results, we reproduce two conflicting experimental measurements, showing 
that these can be understood from a single set of coherent LDA results when thermal effects are 
included. 
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I. INTRODUCTION 



There has been a long-standing controversy about 
the nature of the charge state of the gallium vacancy 
(V r J a ) mediating self-diffusion in GaAs. While exper- 
iments have deduced a preferred charge state ranging 
from -1 and -2 lj to -3 @, Q under N-type or 
intrinsic conditions, most (zero temperature) ab-initio 
calculations favor the —3 ch arge state for Voai see 
Ref. 0SH0H S E3, EH E3, EE O Extracting 
the charge state from experiments is challenging: these 
are usually performed at high temperature, in a regime 
where the formation entropy becomes a significant part 
of the Gibbs free energy. These entropic contributions 
can therefore alter the charge of the dominant species 
as calculated at T = 0, leading to a strong competition 
between the charge states. Moreover, the separate ion- 
ization energy for each state is not directly accessible 
experimentally; only the sum over the ionization free en- 
ergies can be obtained by fitting the experimental data 
over a range of temperatures, using a number of assump- 
tions typically inspired from ab-initio data i or set 
as free parameters in the fit Q]. 

Since the ionization free energies are key for fitting ex- 
perimental data, it is important to evaluate them from 
ab-initio calculations. Because of the high computational 
cost associated with computing vibrational spectra [Tlj . 
however, various assumptions are made typically to sim- 
plify the problem. For example, it is usually assumed 
that the free energy of ionization is independent of tem- 



perature and that ionization levels are constant with re- 
spect to the valence band edge so that the formation 
entropy of all charge defects is taken as equal to that of 
VQ a [Tg] ■ Often, it is also supposed that ionization en- 
ergies follow the thermal behavior of band gap (or the 
conduction band minimum) |17| . which results typically 
into shallower ionization levels with increasing temper- 
ature (within this approximation, for example, heavily 
charged defects in Si become favored at high T [Tslp. 

In this article, we revisit the question of the dominant 
charge state of Vgo, in GaAs, calculating directly, using 
LDA, the ionization Gibbs free energies (F ql / q2 ) for all 
possible char ge s tate transitions, using an approach ap- 
plied to SiCll^ and Si [Hj. We show that, in spite 
of well known LDA limitations, apparently conflicting 
data coming from two experiments performed under dif- 
ferent conditions Q, Q are well fitted. This satisfactory 
fit is obtained for both experimental data by using the 
same set of inputs like the concentration model, the Fermi 
level effect and the computed temperature dependence of 

This article is organized as follows: calculation details 
are presented in section [n] next we analyze the source of 
conflict that arose between different experiments in sec- 
tion lllll presenting at the same time the common concen- 
tration and diffusion models used. Section llVl contains a 
discussion of phonon calculation results followed by our 
fitting procedure and its application on two experimental 
data. The main conclusions are presented in section Ivl 



II. COMPUTATIONAL DETAILS 
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From the computational point of view, the stability of 
V(% a is determined by comparing total energies of super- 
cells with and without defects. In a previous ground state 
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ab-initio calculation using SIESTA, a local-basis LDA 
code 0, we computed the ionization levels of gallium 
vacancies for 0, —1, —2, —3 charges. 

Here, we include the thermal contributions by comput- 
ing the Gibbs free energies in the harmonic limit. Note 
that entropy changes in bulk Ga and As metals is ac- 
counted for (see Refs. [TBI l2p|| as well as Ref. Q for more 
details): 



G = E t , 
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TS vib +p.V + J2»iNi (1) 



where, at constant number of particles and volume, the 
vibrational energy and entropy are given, respectively, by 



ionization levels is an important factor that suffices to ex- 
plain the conflict. Some calculation details are skipped 
in the manuscript, we refer the reader to Refs. 0,0 for 
further details on the model used. 



Concentration model 



Both fitting procedures used by Gebauer et al. Q and 
Bracht et al. [I| are based on the concentration model 
used previously by Tan and coworkers Q. The Ga self- 
diffusion coefficient Dc a is given by the sum of the trans- 
port coefficients of vacancies in various charge states q: 
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and 



where, Co represents the Ga atom density in GaAs 
(C = 2.215 x 10 22 cm _1 ), f q is the diffusion correlation 
3N f^. fy^. K^. factor that contains information about the microscopic 

S v ib = ks E]{"/ — ^I ex P(l — l n [l~ ex P(l — ^r)]}- J um P mechanism, C^ 9 _ g is the thermal equilibrium con- 
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To obtain the Gibbs free energy, the phonon spectrum 
{u>i) is computed for all stable charge states of V^, a using 
a numerical derivative to construct the dynamical matrix. 

The calculation of the dynamical matrix of cells with 
215 sites is a very time-consuming task. Since many 
atomic shells around the vacancy are distorted it is not 
possible to reduce the number of calculations by symme- 
try operations. We thus use the full dynamical matrix for 
the complete supercell. A mesh cutoff of 120 Ry is used 
and the Brillouin zone of the 215-atom cubic supercell is 
sampled at the T-point. The eigen frequencies are then 
extracted at 146 k-points in the Brillouin zone [H^ . 

Except for the splitting of the optical branches at the 
T-point, which cannot be obtained in the harmonic ap- 
proximation, the total phonon density of states and the 
dispersion curve obtained by this procedure for the GaAs 
crystal are in good agreement with experiment and cal- 
culations pH . 122^ . Since we are interested here in the 
entropy difference and not its absolute value, these dis- 
crepancies should not play an important role. 



III. ANALYSIS OF THE EXPERIMENTAL 
RESULTS 

Before discussing our analysis, we first focus on the 
source of discrepancy that exists in explaining different 
experiments. Although the diffusion model and the ex- 
perimental Fermi level dependence used in different ex- 
periments 0,0 are the same, conclusions about the dom- 
inant charge state of Vc a are contradictory. This con- 
tradiction can therefore be due to experimental errors 
or the additional assumptions about the temperature de- 
pendence of the ionization levels used in the experimental 
fits. The following analysis shows that the choice of the 
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is the diffusion coefficient of the 



vacancy Vq^ for q = 0, 1, 2, 3. 

The thermal equilibrium concentration of Vq^ in GaAs 
is defined as 
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where C* is a pre-exponential factor, G^ q is the Gibbs 

V Ga 

free energy of formation of Vq£, T is temperature and ks 
and is the Boltzmann constant. This equation takes into 
account the influence of Asa over-pressure (-Pas 4 ) , often 
used in experiments, on vacancy concentrations. If the 
arsenic over-pressure increases, it affects the stoichiom- 
etry of gallium arsenide by introducing an excess of As 
atoms. With this excess of As atoms, GaAs sample can 
become As-rich and the concentrations of Ga vacancies 
in GaAs are enhanced. 

Different vacancies can introduce energy levels within 
the energy band gap of GaAs. The occupation of these 
energy states depends on the position of the Fermi level. 
Under extrinsic conditions, when the hole or the electron 
concentration introduced by doping exceeds the intrin- 
sic carrier concentration, the Fermi level deviates from 
its intrinsic position. As a consequence, the ratio of the 
charged to neutral vacancy concentrations is changed. 
Vacancies can introduce acceptor levels, i.e., Vq^ intro- 
duces one level located at E xr -\ above the valence band 

V Ga 

edge, Vq£ introduces one at E v -2 and Vq^ introduces 
one at E v -a above the valence band edge. 

The concentration profiles measured experimentally 
contain contributions from all active vacancy charge 
states that cannot be separated. We take advantage of 
the fact that the Gibbs free energy of formation — and 
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thus the equilibrium concentration Cf, 9 - of neutral 

vacancies do not depend on the position of the extrinsic 
Fermi level. Therefore, we can express the total concen- 
tration as function of the ratio of the charged to neutral 
vacancy concentrations: 
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The ratios are given by: 
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where g v -q is approximated to unity (further explana- 
tions and intermediate steps used in this work can be 
found in Ref. 0). 

The concentration of negatively charged vacancies can 
be deduced if we know the position of the Fermi level Ef 
and the ionization energies E v ~ q . Ionization levels do 
not cause much problems as they arc derived from avail- 
able ab-initio calculation or included in the fit. However, 
the Fermi level depends strongly on intrinsic (n,) and 
extrinsic carrier concentrations (n) by the relation: 
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Intrinsic properties of GaAs as function of temperature 
are taken from empirical data , like intrinsic band gap 
Eg, electron and hole concentrations (rii and pi) and in- 
trinsic Fermi level Ej. The band gap of GaAs is supposed 
to vary with temperature following Varshni's |2{| 
empirical expression for GaAs Eg(T) = 1.519 — 5.0405* 
10~ 4 T 2 /(T + 204) deduced from optical experiments |23j . 
The Fermi energy is the important quantity governing 
the equilibrium concentrations of defects. However, it 
changes if any electrically active impurities or dopants 
are present in the material as well as under the presence 
of charged vacancies. The concentrations of these im- 
purities or dopants (Cnonor, G Acceptor) and vacancy con- 
centration in all charge states (C 6 ') must therefore be 

V Ga 

taken into account. The charge balance equation, requir- 
ing that the free electron and hole concentrations n and 
p (which also depend on the Fermi level) must cancel out 
any net charge resulting from the concentrations of all 
positively and negatively charged defects and impurities: 
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The extrinsic Fermi level, for its part, is obtained 
by solving simultaneously Equations El [HI 13 El and| 
self-consistently using the Newton Raphson method [2f 
Once the Fermi level is known, it is used to determine the 
final Gibbs free energy of formation and the correspond- 
ing equilibrium concentrations of all the defects present 
in the structure. Up to this point, the procedure is suf- 
ficient to reproduce the fit of Gebauer et al. and Bracht 
et al. shown in Fig. 7 of Ref. and Fig. 9 of Ref. Q, 
respectively. 



B. Diffusion model 

Few steps are left to fit the measured diffusion coef- 
ficients of Bracht et al. (see the detailed discussion in 
Section 5.2 of Ref. ). Based on the theoretical work of 
Bockstedte and Schefner ^(|, they assume that the en- 
tropy change associated with the Vq^ migration and the 
correlation factor f q are similar for all charge states, i.e., 
D v q ps Dy Ga and ( f q w / = 1). In this approximation, 



Eq.1| simplifies 
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D 



Ga 



1 r< e 1 

Co ^ 



D 



V G , 



E 



CI 



o C v° a 



(13) 



While the extrinsic diffusion coefficient (Dv Ga ) is un- 
known, it is, however, related to the intrinsic diffusion 
coefficient by the relation: 
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The missing quantity in this equation is DG a (ni). Bracht 
et al. use the Ga diffusion coefficient they measured in 
intrinsic GaAs and find Do a , using Eans 1131 1141 and : 



DGa{rii) = 0.64 exp 
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C. Fits of positron-annihilation measurements 

From positron-annihilation experiments, Gebauer et 
al. 0] deduce the Gibbs free energy of formation and 
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FIG. 1: Ionization levels (solid lines) used by Gebauer et 
al. 0, they are taken to have a negative temperature depen- 
dence and vary like the experimental band gap. The position 
of the extrinsic Fermi level is also shown (dash line). 



FIG. 2: Ionization levels (solid lines ) as obtained by Bracht et 
al. from a fit of the experimental data in the temperature 
range 1032K-1350K and are assumed to remain constant. The 
position of the extrinsic Fermi level is also shown (dash line). 



the dominant charge state of Vc a by measuring directly 
vacancy concentrations. Measures were done at different 
high temperatures (900-1450K) by changing dopants con- 
centrations and the chemical potential (stoichiometry). 

The resulting concentration profiles were fitted using 
the model presented above Q using ab-initio ionization 
levels calculated by Baraff et Schliiter 01 a t OK corre- 
sponding to F /- 1 = 0.13E G , F- 1 /- 2 = 0.35E G and 
F~ 2// ~ 3 = 0.49i?G. These ionization levels are supposed 
to vary with temperature like the band gap following 
Varshni's |25| empirical relation. 

Gebauer et al. fitted the concentration profiles for 
N-GaAs with [Te] = 2x I0 18 cm~ 3 at arsenic over pres- 
sure of Pa S4 = 5.6 atm for temperature ranging between 
900 and I450K. The resulting fit (corresponding to Fig. 7 
of their article) is in reasonable agreement with measured 
concentration profiles especially at high temperature, ft 
shows a negative temperature dependence that cannot 
be reproduced unless the —3 charge state dominates all 
other charge states. 

Fig. n shows the variation of ionization levels used as 
function of temperature, the extrinsic Fermi level corre- 
sponding to the the dopant concentration used in this 
experiment is also shown. We can see that this choice of 
the ionization levels favors the stability of the —3 charge 
state for N-type and intrinsic GaAs over a broad temper- 
ature range. 



D. Fits of interdiffusion measurements 

For their part, Bracht et al. p], measured interdiffusion 
in isotope heterostructure 71 GaAs/ nat GaAs exposed to 
arsenic over pressure of Pa s = 1 atm under intrinsic and 



extrinsic conditions (N-type with [Si] = 3x I0 18 cm~ 3 
and P-type with [Be] = 3x 10 18 cm~ 3 ). They measured 
the temperature dependence of the Ga self-diffusion co- 
efficient (DQ a ) at high temperature ranging between 
1032K and 1350K. 

The simultaneous fit of the three curves (N-type, P- 
type and intrinsic curves in Fig. f of Ref. Q) is compli- 
cated and lengthy, especially that the ionization levels of 
Baraff and Schliiter did not give a satisfactory fit. There- 
fore they include the ionization levels in the fitting pro- 
cedure. The best fit obtained corresponds to Hj = 1 .9 
eV and includes two ionization levels only (0/ — f and 
-1/-2) located at F /- 1 = 0.42 eV and F" 1 /" 2 = 0.6 eV 
from the VBM that must remain unchanged with tem- 
perature. 

Figure|21shows the variation of ionization levels used as 
function of temperature. The Fermi level corresponding 
to intrinsic, N-type and P-type GaAs is also shown. By 
analyzing this figure, one can notice that the —2 charge 
state dominates at low temperature, in intrinsic and N- 
type GaAs, then it is reversed to —1 and as tempera- 
ture increases. In this temperature range there is a com- 
petition between different charge states for intrinsic and 
N-type GaAs with no contribution coming from the —3 
charge state. This choice of ionization level favors the sta- 
bility of less charged vacancies as temperature increases. 
For P-type GaAs, the neutral vacancy dominates self- 
diffusion over the entire temperature range. 

In both experiments the dominant part of the tem- 
perature dependence is due to the change of the chem- 
ical potential (i.e., the Fermi-level effect), treated ex- 
tensively by Tan and coworkers 0|. However, the 
above analysis shows unequivocally that only including 
the experimentally-measured thermal dependence of the 
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FIG. 3: Entropy and enthalpy (inset) differences obtained 
after progressive ionization of Vj a . While AH ql/q2 shows lit- 
tle temperature dependence, the non-vanishing entropy dif- 
ference between charge states lead to a linearly increasing 
entropic term, TAS ql/q2 = T(S q2 - S ql ), in the Gibbs free 
energy. 



Fermi level alone is not sufficient to explain the two con- 
flicting experiments at the same time. Figures ^ an( l EH 
show the crucial role played by the ionization level be- 
haviour for determining the charge state of the vacancy 
and that they are the main source of conflict. 
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FIG. 4: Stability domains of Vq in GaAs as function of 
temperature, indicated by the charge state. Each region de- 
notes the temperature and the doping range where a charge 
state dominates self diffusion. The dashed lines show the 
typical position of the chemical potential as function of tem- 
perature for N, P and intrinsic conditions. At temperatures 
above 1000K the competition between different charges be- 
comes significant. 



IV. RESULTS AND DISCUSSION 

The calculated thermal effect on the enthalpy 
A#(9l/g2) and the entropy differences TAS {ql ' q2 ^ be- 
tween the relevant sequential charge states are shown in 
Fig. [3J For all differences, the ionization enthalpy is in- 
dependent of the temperature. In contrast, since the en- 
tropy difference between these states is finite and mostly 
independent of T, the ionization entropy term increases 
linearly with temperature. Moreover, contrary to what 
could be expected, the neutral vacancy has the highest 
entropy so that the formation entropy of Vo a should de- 
crease with an increase in the number of electrons. Im- 
portant changes with temperature occur in the band gap, 
and, depending on the doping type, in the Fermi level. 
This temperature dependance of the Fermi level affects 
the relative energy of neighboring charge states by the 
same amount. It is this effect that is responsible for the 
interesting effect of negative temperature dependance of 
the Voa concentration discussed in Refs. 

While the change in the charge transition level due to 
entropy term is smaller than the shift in the Fermi level 
— the (—2/— 3) transition is moved by about 0.3 eV — its 
effect is charge dependant, affecting the relation between 
the various levels. At OK, Vca ionization levels in GaAs 



are located in the lower part of the band gap and distant 
by 0.1—0.2 eV from each other. Hence, an increase of 0.3 
eV coupled with the Fermi level effect cause important 
changes in the dominant vacancy charge. 

Similar entropic effects are known to affect the prop- 
erties of defects in semiconductors. For example, the 
entropy term of amplitude similar to ours appears to be 
responsible for the lowering the migration barrier in SiC 
with increasing temperature ^5|. It is also associated 
with the destabilization of a metastable configuration in 
GaAs [2]], and the thermal behaviour of the binding 
free energy of acceptor-oxygen complexes in Si . 

These entropic effects are directly associated with 
changes in the phonon density of states (DOS) as elec- 
trons are added to the system, modifying the bonding 
strength in the region surrounding the defect. In the case 
of Vca, the phonon DOS for the neutral vacancy shows 
the highest deviation from the crystalline value. The 
effective impact of the vacancy on the lattice is equiva- 
lent to the application of a local negative pressure: the 
transverse acoustic (TA) modes move from the crystalline 
value by 3 cm -1 to higher frequencies (mode hardening), 
while the LA, LO, TO modes move to lower frequency 
(mode softening) by 2 cm -1 . As the defect is progres- 
sively charged, the DOS tends to recover the frequencies 
of the crystal DOS: TA modes move to lower frequencies 
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while LA, LO and TO move to higher frequencies lead- 
ing to an overall decrease in the entropy with the charge 
state. 

Structurally, this behaviour is associated with the re- 
laxation of the arsenic atoms surrounding the empty site. 
In an unrclaxcd model, As atoms form a tetrahedron 
around the empty site with an interatomic distance of 
3.95 A. As electron are added to the relaxed system, the 
symmetry is conserved but volume of the vacancy de- 
creases. A detailed structural analysis reveals that the 
addition of the first electron does not affect considerably 
the bond length around Va a - As- As bonds decrease pro- 
gressively to 3.37A for and —1 charges and 3.36A and 
3.35A for —2 and —3 charge states. To compensate, since 
the simulation is performed at fixed volume, the Ga-As 
bonds, for the atoms surrounding the vacancy, stretch out 
with the charge; this effect propagates in the direction of 
the link until the borders of the supercell especially on 
the high symmetry axes. The pressure increases locally 
around the vacancy but it is released farthest away from 
it. The phonon density of states suffers from a strong 
competition between the compression and the expansion 
of bonds, and the electron-phonon interaction. 

Finally, it appears that the DOS is under the effects 
of a small progressive compression as electrons are added 
(TA softening and LA, LO, TO hardening) that clearly 
overcomes the proposed phonon mode softening induced 
by electron-phonon interaction studied earlier for GaAs 
in Ref. (2^. Instead, our data show a brutal decrease 
of the phonon vibrational modes as a neutral vacancy is 
introduced in a perfect crystal and a progressive strength- 
ening as electrons are added to the defect, agreeing with 
the experimental observation [3(| that LO(F) phonon 
peak is shifted toward the low frequency side (mode soft- 
ening) and broadens its width depending on the concen- 
tration of vacancies. 

Results for our calculated formation entropies and en- 
thalpies are summarized in Table [I] Formation entropy 
for the neutral state, the only one available in the lit- 
erature, compare well with the ab-initio values using an- 
harmonic approximation by Bockstedte and Schcffler |l6J 
and a recent fit of experimental data by Gebauer et al. |2( 
(further comparisons with earlier works can be found in 
these two references). The formation enthalpy we calcu- 
lated in a previous work Q is in good agreement with re- 
cent LDA calculation by Zollo et al. [3l| using 64 and 216 
supercells, differences do not exceed 0.2 eV for Vq . Wc 
relaxed Vc a in different charge states using spin-polarized 
LDA calculation but no Jahn- Teller distortion has been 
observed. This confirms that the symmetry conserva- 
tion of Vq u of is not a drawback of LDA but a behaviour 
proper to cation vacancies [32l | . To our knowledge no cal- 
culation using GGA has been performed on this defect 
at this point. 

Combining the results presented above, wc obtain the 
free energies difference between the various charge states: 
F /- 1 = 0.22 + 2.74,k B T, F~ x f~ 2 = 0.42 + 1.42fc B T 



TABLE I: Calculated formation entropies Sf and formation 
enthalpies Hf for V^tT 1 ' -2 ' -3 m GaAs are compared to earlier 
calculation and experimental fits (see the text) 

Neutral -1 -2 -3 



S f /kB 7.23 a , 7.3 6 , 9.6 ± l c 4.46 3.04 0.755 

H f (eV) 3.2 a , 2.8 b ,3.2 ± 0.5 C 3.35 3.78 4.423 



a Present work 
b Ab-initio calculation by Bockstedte and SchefHer |16| 
c Fit of concentration profiles by Gebauer et al. [3] 

and i^- 2 /-3 = o.63 + 2.28fc B T. While the formation 
entropies are known to suffer from finite-size effects, en- 
tropy differences should be much more reliable. Forma- 
tion entropy errors due to the finite-size of the supercell 
are expected to be of similar magnitude for all charge 
states and thus cancel out when only differences are con- 
sidered ^5|. Contrary to earlier assumptions 0], the 
actual temperature dependence of the ionization energies 
favors the formation of less charged point defects at high 
temperature. Figure|Hshows the stability domains within 
the experimental band gap for the dominant charge states 
as function of temperature. 

The transition from a charge state to another takes 
place as soon as the Fermi level (under certain doping 
conditions) crosses the ionization energy line, falling into 
a new stability domain. For P-typc, the 0/— 1 level is 
always higher than the extrinsic Fermi energy suggesting 
that the neutral vacancy remains stable at all tempera- 
tures. For intrinsic and N-type conditions, however, the 
charge state of the vacancy is lowered as the temperature 
increases. The critical temperature associated with these 
changes is located in the range of 520-1500K where Vq 2 , 
VQ a and V(Z a start to compete with V^j 3 . 



A. Uncertainties on the calculated Gibbs free 
energy 

The formation enthalpies and transition levels calcu- 
lated with LDA involve uncertainties that can exceed the 
vibrational effects, but that are more likely to cancel out. 
Thus, one must be aware of possible source of errors af- 
fecting the Gibbs free energy of formation. Formation 
entropies suffers from errors associated with the calcula- 
tion of vibrational density of states [3i| ■ For example, 
the microscopic structure of the defect itself, which is 
possibly not obtained accurately. This may be a result 
of the defect-defect interaction of the supercells, but also 
errors from the construction of the pscudopotentials or 
the LDA may contribute. In general, these errors only 
affect slightly the total energy calculation but are more 
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significant for vibrational modes. Moreover, uncertainty 
in the lattice constant which causes the defect to exerts 
pressure on its environment and can modify the equilib- 
rium lattice constant of the defective supercell compared 
to the LDA bulk value. Since we are only interested in 
Gibbs free energy differences we may expect a significant 
amount of error cancellation. Recent tests on phonon 
frequencies of silicon vacancy in bulk silicon [34| showed 
that a cancellation of errors does tend to occur because 
free energy differences depend only weakly on the choice 
of lattice constant. 

Uncertainties on formation enthalpy can come mostly 
from two sources. First, the LDA band gap problem, 
since it is well-known that LDA usually underestimates 
the band gap of materials. The simplest way, used here, 
for dealing with the band gap problem [j] is to shift 
the conduction band states up uniformly by the amount 
needed to reproduce the experimental gap [35| (align the 
conduction band maximum (CBM) with the experimen- 
tal value, at 1.52 eV from the valence band maximum 
(VBM)). Since the LDA can produce similarly large er- 
rors in the energies of the deep defect states, it is also im- 
portant to correct these errors. Therefore, defect states 
with predominantly conduction band character experi- 
ence the same upward shift as the the conduction band 
states themselves, while leaving the defect states with 
predominantly valence band character fixed relative to 
the valence band edge. In our case, V ( | a ionization lev- 
els are calculated by total-energy energy differences, the 
correction we use here leaves these ionization levels un- 
changed since they have a valence band character. Many 
other correction schemes can be used with good efficiency 
as discussed in the work of Castleton and coworkers [36| , 
we discuss below the impact of such corrections on the fit- 
ting procedure. With current computational power limi- 
tations, a full GW calculation, which would correct LDA 
errors, is unfortunately not possible for the large super- 
cells needed for defect studies. Therefore, LDA calcula- 
tion is the best choice currently since it provides a good 
estimate of the electronic ionisation levels that enables 
us to fit efficiently experimental data. Second, supercell 
corrections could also introduce additional errors on the 
calculated formation enthalpy. In this work, total ener- 
gies of charged defect are corrected using the first order 
Makov-Payne method H^. Many articles [M III El 
have recently shown that Makov and Payne correction 
can be overestimated for highly charged defects. 



B. Fitting experimental measurements 

We revisit the analysis of two apparently conflicting 
results by looking at diffusion data of Bracht et al. and 
concentration profiles of Gebauer et al. We first repro- 
duce exactly the fits of Gebauer et al. and Bracht et al. 
using exactly the same conditions described in the respec- 
tive original papers (dash lines in Fig. [S] and squares in 
Fig. inj an d explained above. We then use the Gibbs free 
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FIG. 5: Experimental measurement (filled circle), and fit 
(dashed line) for temperature dependence of the vacancy con- 
centration, as reported by Gebauer et al. . Our calculation 
including —2 and —3 contributions (solid line) is shown to- 
gether as well as separately (dotted line). 



energy of formation and all resulting ionisation levels cal- 
culated by LDA to generate our fits. The contributions 
coming from each charge state is shown and the total 
concentration and diffusion coefficient are plotted with 
solid lines in Figs. [S] and [fl] 

Using our calculated values, with only one free param- 
eter to adjust the prefactor (C* in equation |5J) , and —2 
and —3 charges, we reproduce the experimental data of 
Gebauer et al. with a good accuracy (Fig. |3J) showing 
that the —2 charge state is the main contributor to self- 
diffusion at high T. 

Self-diffusion coefficient measured by of Bracht et al. 
are also efficiently fitted. Comparing our calculated Ga 
self-diffusion parameters with experimental data and the 
corresponding fits 0], with an adjustable prefactor C* we 
notice (see Fig. EJb)) a strong competition between —2, 
— 1 and charge states for intrinsic GaAs while —3 does 
not contribute to the total concentration. Moreover, a 
progressive charge reversibility occurs: Vq^ dominates 

at high temperature until T^ 2/ ~ 1] = 1230K, then 

dominates until Tj = 1345K. Even if details about 
the location of transition temperatures are different, the 
total diffusion coefficient is very well reproduced with 
our model too. Contrary to what was proposed by Tan 
et al. Q, our results support the experiment, indicating 
that under intrinsic conditions V^T 3 plays very little role. 

For Si-doped GaAs, while Bracht et al. find that the 
concentration decreases with temperature for V^T a 2 but 



shows no contribution from the —3 charge states, we see 
a more complex situation (Fig.|HJa)): self-diffusion is me- 
diatcd by V G * at temperature below ' ' = 1150K, 
while above this transition temperature —2 charge state 



is dominant until T 



(-2/-1) 



N 



1450K. Finally, for P-type 
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FIG. 6: Gallium self diffusion constant versus inverse tem- 
perature for (a) N-type (b) intrinsic (c) P-type. Theoretical 
results for D q Ga (dashed lines) show the contribution of each 
charge state. The transition temperature lie at the crossing of 
the curves. The total Dca (solid) is in good agreement with 
the measured interdiffusion data (square) and fit by Bracht 
et al. (circle). 



GaAs (Fig. Etc)), we obtain results that are very simi- 
lar to Bracht: the neutral vacancy is dominant over the 
whole temperature range. 

As is seen in Figs. [S] and [ft] our model reproduces very 
well all data at high temperature for different doping con- 
dition. The agreement is weaker for N-typc GaAs below 
HOOK, however. In this regime both experimental data 
sets lie below our predictions. This is likely due to the for- 
mation of vacancy complexes and slow Ga-self-diffusion 
which appear to delay the establishement of equilibrium 
condition, as noted by Gebauer et al. Q. Thus, the mea- 
sured concentration is a lower limit of the calculated equi- 
librium concentration. 

A gap or supercell correction that applies the same 



rigid shift to the ionization levels would not change the 
observed sequential charge reversibility, caused by en- 
tropic effects, but only the location of transition tem- 
perature. In this particular case, if the ionization levels 
are shifted upwards, the contribution of —2 charge state 
increases. Thus, the remaining discrepancies that ap- 
pears between experiment and theory visible in Fig.EJa) 
for the low temperature limit can be significantly atten- 
uated. At a first insight, we propose to use the simplest 
LDA corrections, if entropy effects alone give unsatisfac- 
tory results, one can go beyond by using more advanced 
gap correction schemes. 

V. CONCLUSION 



The present work shows that the deformation potential 
and the electron-phonon interaction induced by charged 
vacancies affect the phonon spectrum in a much more 
complicated way than was thought previously [Til Il7j . 
Consequently, the picture about ionization levels ob- 
tained at OK Q is not applicable in the experimental 
range as ionization free energies are temperature depen- 
dent. We have found that the entropy of formation cal- 
culated with LDA contributes to lower the free energy 
of formation at high temperature when associated with 
empirical thermal effects on the Fermi level. Thus, the 
formation of less charged Ga vacancies mediating self- 
diffusion in GaAs becomes much more favored at high 
temperature. In spite of LDA limitations, the resulting 
coherent Gibbs free energies of formation are a possible 
solution to the apparently contradictory experimental re- 
sults about the dominant charge state for vacancy diffu- 
sion in GaAs. In addition to resolving this controversy, 
our results underline the need to include thermal effects 
on the Fermi level together with ionization entropies and 
ionization enthalpies. Even though the Fermi level effect 
dominates the correction term, the smaller contribution 
coming from ionization entropy is crucial for a correct 
interpretation of experimental data. 
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